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Abstract. We propose a fully nonlinear framework to construct consistency relations for 
testing generic cosmological scenarios using the evolution of large scale structure. It is based 
on the covariant approach in combination with a frame that is purely given by the metric, the 
normal frame. As an example, we apply this framework to the ACDM model, by extending 
the usual first order conditions on the metric potentials to second order, where the two 
potentials start to differ from each other. We argue that working in the normal frame 
is not only a practical choice but also helps with the physical interpretation of nonlinear 
dynamics. In this frame, effective pressures and anisotropic stresses appear at second order in 
perturbation theory, even for "pressureless" dust. We quantify their effect and compare them, 
for illustration, to the pressure of a generic clustering dark energy fluid and the anisotropic 
stress in the DGP model. Besides, we also discuss the effect of a mismatch of the potentials 
on the determination of galaxy bias. 
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1 Introduction 

Current cosmological observations provide evidence for a recent onset of a phase of accel- 
erated expansion of the Universe [1-5]. Under the assumption of large-scale homogeneity 
and isotropy, the observations can only be described in the framework of General Relativity 
(GR) by invoking some form of Dark Energy (DE) with negative pressure. Apart from the 
Standard Model (SM) of particle physics, the preferred model of the content of the Universe 
also assumes the presence of Cold Dark Matter (CDM) and DE in the form of a cosmolo- 
gical constant (A) and is therefore named ACDM. The dynamics of the Universe is based 
on the Einstein equations and on the assumption that the average geometry is described 



- 1 - 



by a Friedmann-Lemaitre- Robertson- Walker (FLRW) metric. The formation of structures is 
then modelled using cosmological perturbation theory [6-11] which describes the evolution 
of small fluctuations in the energy-momentum content and the metric. 

The ACDM model comes with a minimal set of parameters that are increasingly well 
constrained by the data and, to date, no significant departures from its predictions on cosmo- 
logical scales have been found [3, 12, 13]. However, the difficulty in explaining the required 
tiny value of A has motivated many alternative scenarios with dynamical DE fields or modi- 
fications to GR aka Modified Gravity (MG) [14-24]. Given the existing tight constraints, it 
is essential to understand very accurately what the theoretical predictions are and how to 
interpret the observables, in order to maximise the benefits of the future data. A possible 
falsification of ACDM would then be the first step towards a complete understanding of the 
dark energy. In practice, we can choose between comparing ACDM to a large representative 
number of alternative scenarios, parameterise departures from it in a sufficiently flexible (but 
still economical) way [25-39], or identify key relations whose breakdown would hint towards 
a falsification of the standard paradigm [28, 40, 41]. All three strategies are useful in their 
own right: the third option is initially most powerful for falsifying or confirming ACDM, 
while the second can be helpful as a guidance for the first, more traditional option. In any 
case, exquisite precision in the model predictions and understanding of the observables are 
needed to make decisive statements. 

In this work we propose and start to explore a program for constructing fully nonlinear 
consistency relations between geometry and matter content, not assuming GR but allow- 
ing for generic modifications of the field equations. This can be seen as a (nonlinear) step 
beyond the popular modified growth parameterisations [25-39]. Such relations can either 
be implemented as consistency checks to falsify ACDM (or some other paradigm), or can 
be parameterised to explore a possible departure from ACDM. To construct the consistency 
relations we make use of the fact that at late times matter (baryonic and dark) is well de- 
scribed by a pressureless perfect fluid (dust) on cosmological scales. The absence of isotropic 
and anisotropic pressures in the matter rest frame directly implies two consistency relations 
between the geometry and the physics that gives rise to the late time accelerated expansion. 
However, we point out some important subtleties that arise due to the choice of frame and 
that can become relevant on nonlinear scales. 

We first formulate the consistency relations in the exact nonlinear theory employing 
the covariant approach [42-47]. Then, we discuss how the conditions come about after 
matter-radiation equality in first and second order cosmological perturbation theory, using 
the longitudinal (or conformal Newtonian) gauge. At first order, the resulting conditions 
are well known. They simply tell us that the two metric potentials are equal, <f>\ = rpi, 
and their evolution is governed by a simple second order differential equation that can be 
solved exactly: the Bardeen equation [6]. At second order, the total measured pressure 
and anisotropic stress depend on the observer frame that is used to define these quantities, 
and this is where the subtleties come into play [48]. The pressure measured in the frame 
comoving with the matter is zero by definition (in analogy to the Lagrangian picture in the 
Newtonian approximation, see for instance [49, 50]). However, in a non-comoving frame, 
effective pressures and anisotropic stresses are induced by the relative velocity field, even if 
the source is a pressureless perfect fluid [48]. As we shall discuss in detail, there are good 
reasons to work in a specific non-comoving frame, called the normal frame. This frame is 
orthogonal to the surfaces of constant time and is given purely by the geometry, i.e. the 
(perturbed) metric. An added advantage of this framework is that the density and pressure 
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of matter (and of any other component) in the normal frame are then immediately related 
to the geometrical pressures that can be read off from the Einstein tensor. This means that 
projecting on the normal frame lets us cleanly separate the geometrical fluctuations from 
those of the stress-energy content. Consequently, the consistency relations in the normal 
frame take a particularly simple form. The change of frame from comoving to non-comoving 
induces effective matter pressures that we quantify. We find them to be small in comparison 
to typical non-standard sources of pressure and anisotropic stress in models beyond ACDM. 

The outline of the paper is as follows. First, in section 2, we review the effects of 
changing the observer frame from comoving to non-comoving, and discuss why the normal 
frame is particularly useful. In section 3 we state the consistency relations in the exact 
nonlinear case, discuss the connection between the covariant and the metric perturbations 
approach and review the first order consistency relations. In section 4 we derive the second 
order consistency relations and compute the effective matter pressure and anisotropic stress. 
Section 5 is devoted to the discussion of a number of typical cases where the effective matter 
stresses could potentially lead to wrong physical interpretations of the data. We consider the 
cases of a general clustering DE fluid and the DGP (Dvali-Gabadadze-Porrati [17]) model 
as a prototypical example of MG, because for both their respective intrinsic pressure and/or 
anisotropic stress perturbations could, in principle, be confused with the contributions from 
matter. Further, we discuss the impact on the determination of the galaxy bias from weak 
lensing and galaxy clustering. We summarise our work and conclude in section 6. We include 
three appendices. In appendix A we state our conventions on notation and units. Appendix 
B briefly reviews the covariant approach and appendix C gives more details on the metric 
perturbation calculations. 

2 Covariant approach and change of frame 

The covariant approach is very useful to formally describe the full nonlinear dynamics in terms 
of physical variables that are covariantly defined as projections with respect to a family of 
observers, i.e. a frame. In this section, we first discuss the notion of observer frame in the 
covariant approach and its relation to the choice of coordinates and gauges in cosmological 
perturbation theory. Then, we focus on the effect of changing the frame in the case of pure 
pressureless dust, which forms the basis for stating the fully nonlinear consistency relations 
in the next section. 

2.1 Frames and gauges 

The covariant approach is based on the choice of a unit time-like 4-velocity that defines a 
1+3 splitting of space-time and a unique decomposition of tensorial quantities into irreducible 
scalars, and projected vectors and tensors. Such a velocity field defines an observer frame 
and is generically denoted by (with u a u a = — 1). The general decomposition of an energy- 
momentum tensor 

TV = pu"u" + pV^ + 2q^u v ^ + tt^ (2.1) 

yields the dynamical quantities (or fluid variables): energy density p, pressure p, energy flux 
and anisotropic stress Tr^ as measured by a family of observers with 4-velocity u^. The 
tensor V^ v = g^ v + u fl u u is a projector on hypersurfaces orthogonal to u M and we follow the 
convention that indices in round brackets are symmetrised. For more details on the covariant 
formalism and the relevant definitions we refer the reader to appendix B.l. 
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Let us remark that an observer frame is not a choice of gauge (or coordinates). The 
dynamical quantities in a given frame can in general be described in any gauge. However, 
the coordinates in which a specific frame takes its simplest form define a gauge which is 
intrinsically related to the choice of frame. For example, the comoving frame of matter 
(defined by the matter 4- velocity) takes its trivial form ujj, = a~ 1 5^ in the synchronous 
comoving gauge where the matter peculiar velocity vanishes. In this sense, a choice of frame 
prefers (but does not demand) the specific gauge in which the frame is trivial. 

Choosing the comoving frame to do perturbation theory in curved space is analogous 
to using the Lagrangian picture in the Newtonian approximation, where the observer moves 
with the fluid flow. On the other hand, non-comoving frames can be seen as the analogues 
to the Eulerian framework in the Newtonian limit, where the observer and the coordinates 
are independent of the fluid and the properties of the fluid change at fixed coordinates in 
3-space as time flows. However, in curved space-time there are infinitely many frames that 
are not comoving with, say, the matter. We can always translate one choice of frame into 
another but, for practical purposes, the question is how to choose a useful one. A hint can 
be found from the fact that non-comoving frames that see neither shear nor vorticity, called 
quasi-Newtonian frames, are those that provide the most direct connection between general 
relativistic and Newtonian cosmology, as argued e.g. in [51, 52]. 

For any given metric, the lapse and shift, N and N l , of the Arnowitt-Deser-Misner 
(ADM) 1+3 formalism can be used to define a covariant vector field normal to the spatial 
hypersurfaces, = —N5®, = A^ _1 (l, —N l ) [53]. In a perturbed FLRW metric this 
simply becomes 

- ^ (2.2) 

where n a n a = — 1 and rj is the conformal time. We will refer to the vector as the normal 
frame. Clearly, this frame is completely defined by the metric. The coordinates in which the 
normal frame has vanishing spatial components, ri 1 = = rii, are those of the longitudinal 
gauge. Moreover, the shear and the vorticity seen by vanish in this gauge, which means 
that the normal frame is quasi-Newtonian. The zero-shear condition for can in fact be 
used to define the longitudinal gauge at first as well as at second order in perturbations, see 
e.g. [54]. In the appendix C we provide details on the properties of the normal frame and its 
kinematic quantities in the metric perturbation approach. 

The normal frame is a non-comoving frame that is useful as a computational and 
conceptual tool for the following reasons. The evolution equations for matter in second order 
relativistic perturbation theory in the normal frame coincide with the Newtonian hydro- 
dynamic equations in the Eulerian frame, see e.g. [49]. Furthermore, the longitudinal gauge 
variables match the output of N-body simulations, except for the density perturbation that 
needs a gauge-correction on large (linear) scales where the Newtonian approximation starts 
to break down [55, 56]. Finally, we will show in section 4.2 that the contributions of the 
stress-energy to the source terms of the evolution equations of the second order metric po- 
tentials are directly related to the pressure and anisotropic stress in the normal frame. This 
provides a nice physical interpretation of the normal frame quantities. 

2.2 Change of frame 

Fluid variables computed in a given observer frame can be translated into those seen in 
another frame fully nonlinearly, see e.g. [46, 47]. The transformations are given by the 
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relative velocity between the two frames, v^. For the fluid variables all transformations are 
second order in the relative velocity, except for the energy flux. 

Here we focus on the transformation from the (comoving) rest frame of matter, denoted 
by Um, to the (non-comoving) normal frame n^. The general case can be found in appendix 
B.2. By definition, the rest frame (also called energy frame) is the one in which the energy- 
flux vanishes [44]. We denote the fluid variables measured in the rest frame with an asterisk 
and therefore qm = 0. The rest frame velocity decomposes with respect to the normal frame 
as 



< = lm K + O (2-3) 

where 7 m = (1— ^) _1/ ' 2 (with = D ma ti") is the local Lorentz factor of the transformation, 
and the relative velocity Vm is orthogonal to n M such that n Q w^ = 0. The fluid variables as 
measured by the normal frame are related to the rest frame variables via 

Pm = l 2 m p*m ( 2 - 4 ) 

Vm — ^lm v rnPm ~ ^ v mPm 

Qm = Pm v m + ImVmPmVrn = PmV^ (2-6) 
n m ~ lmPm u m u m ~ PmV m u m \*"*) 



where we used the transformation of the energy density (2.4) to express the effective pressure, 
energy flux and anisotropic stress in terms of the energy density in the normal frame and the 
relative velocity. The non- vanishing stresses present in the normal frame, p m and TTm , are 
therefore related to the motion of the matter particles. It is also instructive to realise that 
the pressure corresponds to the difference in energy density measured in the two frames 

Pm= o (Pm ~ P* m ) , ?C = (pm ~ P* m ) ( 2 -8) 

^ V m 

and Pm — Pm = v mPm- In other words: the trace of the energy- momentum tensor, = p — 3p 
is invariant under the change of frame. 

To get some more intuition, let us see how the non-comoving density, p m , differs from 
the comoving density p* m . The mean pressure is related to the relative difference of the mean 
densities, defined as 

v = (Pm(x)) ~ (p* m (x)) = (3p m (x)) 
(Pm{x)) (Pm(x)) 

In second order perturbation theory, this turns out to be the dispersion of the relative velocity, 
see section 4.1. In terms of the two-point function £(r) = (p m (x) p m (x + r))/(p m (x)) 2 — 1 we 
find that the difference between the comoving and the non-comoving correlation function is 
given by 

r (r) - f(r) = (2u - , 2 ) (t(r) + 1) + {3pm{x) 3Pm(X + 0) 7 2{ ' m{x) * Pm(X + r)) (2.10) 

which nicely shows that this difference is mainly due to the velocity dispersion with correc- 
tions involving the pressure auto-correlation and pressure-density cross-correlation. Without 
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further investigations we cannot say which is the precise correlation function actually ob- 
served. But since we do observe effects from peculiar velocities, like redshift space distor- 
tions, it appears unlikely that observations are not sensitive to the effective pressure and 
anisotropic stress. 

To summarise, the fluid variables measured in the two types of frames coincide at first 
order in perturbations (except for the energy flux that is linked to the relative velocity). 
However, at second and higher order there is an intrinsic difference as effective stresses arise 
in non-comoving frames. We will quantify p m and Hm in second order perturbation theory 
in section 4.1. As we will see in the next section, the presence of the matter pressure in non- 
comoving frames translates into new terms in the nonlinear consistency relations between 
geometry and matter. For further insight into the role of matter pressure in non-comoving 
frames, see also [48]. 

3 Exact consistency relations and linear theory 

The currently most promising approach to learn about the physical nature of dark energy 
requires not just an accurate determination of its equation of state, but also of its pressure 
perturbation and anisotropic stress. In a ACDM universe, the pressure and the anisotropic 
stress vanish in the matter rest frame (apart from a small contribution from relativistic 
particles that is negligible at late times). Any deviation from zero would therefore be nor- 
mally interpreted as a sign of non trivial dark energy dynamics. Let us now see how this 
consistency check comes about in the covariant approach and how its precise statement ac- 
tually depends on the observer frame. Then we will make the connection to the metric 
perturbations approach and discuss the consistency relations in linear theory. 

3.1 Fully nonlinear ACDM consistency 

We start with a minimal set of assumptions before restricting the discussion to the ACDM 
scenario. We assume that matter follows the geodesies of a metric that gives rise to an 
Einstein tensor G^ v '. Matter refers here to baryonic matter, cold dark matter, radiation and 
neutrinos; and it is described by the total (matter) energy- momentum tensor T^ v '. Then, 



is interpreted as the effective energy-momentum tensor of the dark physics causing the late 
time accelerated expansion. X^ v can represent a cosmological constant, any additional dark 
energy fields, modifications to GR, variations of Newton's constant or effective back-reaction 
effects. We refer to X^ v collectively as the dark energy. The only assumption so far is 
that there is (effectively) a single metric that describes space-time and gravity. If we also 
assume that the theory is free of torsion, the Bianchi identities of the Einstein tensor read 
V ' a G a ^ = and therefore yield V Q (T a/i + X a ^) = 0. Only under the further assumption that 
T^ u is covariantly conserved on its own we would also be allowed to write that V a X a ' i = 0, 
but we do not need to make this restriction here. 

The path to get the full nonlinear consistency equations is to decompose (3.1) just like 
the energy- momentum tensor in (2.1). In order to learn about the dark energy, we could 
then project it using its own 4-velocity which is, though, unknown. Conversely, if we were 
to read off the dark energy properties in the matter rest frame, we would decompose (3.1) 
with respect to Um and deduce 



= M 2 pG^ v - 



V 



(3.1) 




7T 



X = 7r G 



(3.2) 



- 6 - 



However, the geometric pressure and anisotropic stress measured in the matter rest frame, 
Pq and t^*q U ■, are nonlinear combinations of matter and geometry variables that are hard 
to disentangle. Projecting on the normal frame will separate geometry and matter neatly. 
Therefore, we propose to decompose (3.1) in the normal frame and apply the transform- 
ations (2.5) and (2.7) to find 

Px = PG ~ \v 2 m p m , 7T^ = - pmV^V^ . (3.3) 

Since the normal vector is defined exclusively by the metric, the geometric pressure and 
anisotropic stress measured in the normal frame, pc and -n^f , are purely given in terms of 
the geometry, i.e. the metric potentials and the scale factor (for FLRW). Therefore, these 
quantities can in principle be reconstructed from purely geometrical observations such as 
Weak Lensing (WL), the Integrated Sachs- Wolfe (ISW) effect, the motion of test 'particles' 
like galaxies, and the background expansion measurements. The contributions to the dark 
energy variables from the matter pressure and anisotropic stress can be reconstructed from 
observations of clustering, peculiar velocities and growth of structure. The exact nonlinear 
ACDM consistency relations are now simply deduced by requiring px = —MpA and ir^ = 0: 

- MpA = ( PG - \v 2 mPm ) , = ( tv g u - p m v^v v S] (3.4) 

V 6 J ACDM V / ACDM 

Any breakdown of these relations would falsify the ACDM paradigm, and the dark pressure 
and anisotropic stress deduced from (3.3) would give information of the physics beyond 
ACDM. This is one of the main results of this paper: the velocity contribution to these 
consistency relations is absent in linear perturbation theory. The commonly used ACDM 
consistency relations at first order are therefore broken by this contribution. 

3.2 Comparison with modified growth approaches 

Let us now briefly connect this framework to the modified growth parameterisations of [25- 
39]. These are sometimes also referred to as Parameterised Post-Friedmannian approaches. 
First of all, notice that these formalisms are usually restricted to linear scalar perturba- 
tions, while the relations (3.3) and the ACDM consistency conditions (3.4) are fully general, 
nonlinear and are not restricted to scalar degrees of freedom. 

The ratio between the gauge invariant scalar metric potentials (ft an d ip is used as an 
indicator of the anisotropic stress of X^ u ', see section 5.2. In addition, the contribution from 
X^ v to the Poisson equation is absorbed into a time and scale dependent modification of 
Newton's constant, defined as G e s = G(l + Apx / Ap m ), where Apx and Ap m are the density 
perturbations in the comoving orthogonal gauges of matter and dark energy, respectively. The 
nonlinear counterpart to this approach in the framework that we propose is straightforward: 
the anisotropic stress is measured through the second equation of (3.3) and the modified 
Newton's constant is given by G g = G(l + Sp x /Sp m ), where now 5p x and 5p m are the 
fluctuations measured in the normal frame and we can write 

Gcflf/G = 5p G /5p m . (3.5) 

G e ff is usually defined in Fourier space. Notice however that when taking G e g to be a ratio 
between perturbed quantities it makes a difference if it is defined in real or Fourier space. 
On the other hand, our consistency relations are defined in real space and can directly be 
translated into Fourier space. 
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3.3 Metric perturbations and the normal frame 

In this section we provide the connection between the covariant approach and second order 
metric perturbations. We choose to work in the generalised longitudinal gauge, where the 
spatial components of the normal frame vanish and there is zero shear. For a description of 
the second order longitudinal gauge for scalar quantities see for instance [54] . 
The homogeneous and isotropic background evolution is characterised by 

PG = 3Mj>a- 2 n 2 = - Pm + Px (3.6) 
PC = -M 2 P a- 2 (m' + U 2 ) = p m + px (3.7) 

where the prime is the conformal time derivative and the conformal Hubble parameter is 
% = a! I a. We can define the geometrical or total equation of state parameter as 



PG 3 V K 



fo = ^=-i(l + 2;£r • (3-8) 



Measurements of the background evolution via standard candles and rulers only give T~L(a) 
and not p m . As a consequence, Pg = Pm + Px and pc = Pm + Px cannot be disentangled 
without further assumptions [57]. 

Let us focus on scalar perturbations only, even at second order 2 . The shear seen by the 
normal frame vanishes when the off-diagonal parts of the metric are set to zero by choosing 
the gauge functions appropriately. At first and second order we are only left with the usual 
scalar potentials in the diagonal of the metric, tp = tpi + hip2 an d <p = (pi + ^>4>2, where the 
subscripts indicate the order of the perturbation. These metric perturbations coincide with 
the Bardeen potentials in the longitudinal gauge. For practical reasons that will be clear 
later on, we choose to work with the variables 

fn = ^ (<t>n + VVi) , U n = 4> n -1p n (3.9) 

where again the subscript n denotes the order of the perturbation. In the following, wherever 
this label is omitted we mean the first plus the second order, according to the following 
convention: tp = pi + \p>2- The perturbed metric then reads 

g^dx^dx" = a 2 {- (1 + 2ip - II) dry 2 + (1 - 2cp - II) 5 i:j dx l dx j } . (3.10) 

The so called Weyl potential, ip, is related to the WL observables and the ISW effect in 
the Cosmic Microwave Background (CMB). The potential II is related to the anisotropic 
stress of the theory and quantifies the mismatch between the original potentials ip and cj). 
This potential, that is often called gravitational slip, is not directly related on its own to 
specific observables, although it can in principle be reconstructed from combining WL ob- 
servables with galaxy clustering, as we discuss briefly in section 5.3, or with peculiar velocity 
measurements, e.g. from redshift-space distortions; see also [26, 28, 41]. 



1 In ref. [54] the normal vector and the kinematic quantities are given without specifying the gauge. Note 
that in our notation the roles of <j> and ip are swapped with respect to theirs. 

2 Clearly, this treatment is not complete as second order vectors and tensors are sourced by products of 
first order scalars and tensors. Since the scalar sector dominates, we assume, for the scope of this work, that 
the effect of vectors and tensors is negligible. 
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In these metric variables, the normal vector on constant conformal time hypersurfaces 

reads 



2 



nf = a-^|i- v + In + |( w -in x j I (3.11) 

n» = -a§ |l + ^-in-i (vi-^ni) I (3.12) 

and the acceleration and isotropic expansion experienced by can be found in appendix 
C.l. The shear and the vorticity vanish in this gauge by construction. The relative velocity 
between the normal frame and the matter rest frame is defined in (2.3) and thus satisfies the 

1 /2 

equation Vm = (l — i^J n y — n M . This can be solved perturbatively to find 

v% l = a- 1 6?Vi, v m ^ = aSl l [V mi -(2cp 1 + U 1 )V mU ] (3.13) 
where V m i = = V^ nl + iV^ 2 is the usual velocity perturbation of the matter 4- velocity 



u 



ii 



a _1 (l^, KJi); which satisfies = —1. Finally, the perturbed fluid variables are 
written in the form p = p + 5pi + \5p 2 (and analogously for the pressure, energy flux and 
anisotropic stress perturbations). 

Notice that the definition of the perturbation variables depends on the observer frame. 
Therefore, at second order our perturbation variables are non trivially related to the vari- 
ables in the usual formalism where the rest frame fluid variables are chosen. As usual, the 
density and pressure perturbations are not gauge invariant. One may define gauge invariant 
variables by taking the covariant spatial gradients of the density and pressure [44]. The 
reason is that variables that vanish at the background level are gauge independent at the 
lowest non- vanishing perturbative order - a consequence of the Stewart- Walker lemma [58] . 
This also means that 5p m2 is gauge independent because p m = = 5p ml , and analogously 
for the anisotropic stress. In fact, <5p m2 and dir m 2 coincide with particular gauge invariant 
combinations in the longitudinal gauge, see e.g. [48]. 

3.4 Consistency relations at first order 

We now compute the geometrical fluid variables at linear level by decomposing G^ u in the 
normal frame. Focusing on the late time evolution of ACDM, the absence of pressure and 
anisotropic stress perturbations in matter and the form of X^ u = —MpAg^, lead to the 
simple linearised versions of the consistency relations (3.3). These linear relations reflect the 
fact that the two scalar potentials are equal and their evolution is given by the well-known 
Bardeen equation for tpi, which we can solve analytically. 

The geometrical fluid variables in the normal frame for linear perturbations are found 

to be 

u m 2 2 1 m m 
2 / u/, <p" 3 m n" v 2 rii\ 

Spa = -fPc -3 TO + 3 £ + ^ + + ^ + ^ - w ) (3.15, 



5p G1 = -2-p G ( <p x + ^ - - -Hx + ^ - ) (3.14) 



2 5^ 
5q ^ = —Z pG «H di 



f\ Itt n 'l 



(3.16) 



1 h % ^^ v / 1 \ 

S*& = - 3 PG (didj - 3^V 2 ) ni . (3.17) 
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Well after matter-radiation equality, the matter can be modelled as pressureless dust and 
the perturbations in the radiation can be neglected. At first order, we have from (2.5) and 
(2.7) that 5p ml = and 5ir m i = 0, while 5p ml = p m 5 ml and 5q^ = pmV^a' 1 8? , where 5 m \ 
and are the first order density contrast and velocity perturbation, respectively. Notice 
two interesting consequences: first, if the dark energy is not a cosmological constant, its 
isotropic and anisotropic pressure perturbations, 5p X i and Sttxi, can directly be inferred from 
the metric potentials via (3.15) and (3.17), respectively. Conversely, we can also conclude 
that for any metric inferred from observations there exists a pressure perturbation 5p X \ 
and an anisotropic stress Sttxi, given by the equations above, that describe a dark energy 
component which is able to reproduce the observations. Second, if ACDM is the correct 
model, the geometric pressure and anisotropic stress perturbations vanish identically and the 
exact consistency relations (3.4) reduce to 

ACDM IIi = and if" + 3?Vi - 3w G H 2 <pi = . (3.18) 

The resulting equation for ipi is the well-known Bardeen equation for dust after matter- 
radiation equality. Note that no spatial derivatives are present, which means that ipi evolves 
in a scale-independent way. After matter-radiation equality, we can write the geometric equa- 
tion of state in ACDM as wq = — ^a/ (^a + ^m a_3 )> where f2« is the present energy fraction 
in the z-component. Then, the equation for </?i can be solved analytically, see appendix C.2. 
Neglecting the decaying mode, which soon becomes irrelevant, we find in Fourier space 

<Pl = C k T{-a 3 ) (3.19) 

where a = (f2A/f2 m ) 1//3 a is the scale factor normalised to 1 at A-matter equality, and 
J-{y) = 2-Fi(l/3, 1; 11/6; y) is the ordinary hypergeometric function. J-(— a 3 ) evolves from 
1 in matter domination to about 0.76 today and describes the suppression of the growth 
of perturbations due to the cosmological constant. The integration constant is fixed by 
the initial conditions that are set well after radiation-matter equality. It can be estimated 
analytically by solving the perturbation equations of a radiation-matter mixture from the 
end of inflation until matter domination, neglecting the anisotropic stress of radiation, see 
e.g. equations (2.216) and (3.65) of ref. [59] 3 

3 / — ~ — ~ 7T f 1 for k < k s , 

C^ = --J2Ti 2 k-^VAk)x { _ 3.20 

5V cv > \2>{k s /k) 2 hi{k/k s ) for k^k s K ' 

where V^(k) = A s (k/k p ) ris ~ 1 is the primordial dimensionless power spectrum of the comoving 
curvature perturbation (, with A s ~ 2 x 10 -9 and n s ~ 0.96 at k p = 0.05/Mpc. We have 
introduced the sound horizon at matter-radiation equality, k s = ^/3/r] cq , with conformal 
time at equality, r/ eq = 2(2 — \/2)/H eq , and H cq = HQ^} m ^2/^} r « 0.010 /i/Mpc, such that 
k s ss 0.015 h/Mpc. 

We can derive the matter density contrast and the velocity perturbation from <p% by 
means of (3.14) and (3.16), see appendix C.2 for details. In the next section we will use the 
velocity perturbation to compute the effective matter pressure and anisotropic stress. We 
find 

Vi 1>k = -iV(a)^C k (3.21) 
Ho 



3 The coefficient for large k is actually given as 20 in ref. [59], but we find numerically that the correct 
coefficient is 3. 



- 10 - 



which is valid on all scales. Its time evolution is given by 



V(o) 



2 


"a(l + a 3 )~ 


1/2 


3 







^(-a 3 ) - 3alF( 



-a 



(3.22) 



where T(y) = (d/ dy)F(y). 

Let us notice that the effortless way in which the Bardeen equation was derived here is 
due to the application of the consistency relations of the geometric pressure that we intro- 
duced in the last section. Next we will see how this is generalised at second order and how 
the effective matter pressures come into play. 



4 Second order perturbations 

In this section, we derive the consistency relations in second order perturbation theory. When 
going beyond linear perturbation theory, the choice of the observer frame becomes important. 
As discussed in section 2.2, the transformation from the rest frame to the normal frame gives 
rise to an effective pressure and anisotropic stress proportional to the square of the relative 
velocity, even for a pressureless perfect fluid. We compute and quantify these contributions 
in ACDM such that they can be properly taken into account when applying the second order 
consistency relations. 

4.1 Effective pressures 

The effective pressure and anisotropic stress induced by the change of frame are given in 
exact form in (2.5) and (2.7). We use the perturbed relative velocity from (3.13) to find 

2 _ 

SPm2 = ^PrnV^VmU (4.1) 

6*% = 2p m jCi^i - l^V'aVmik} a- 2 S?5» . (4.2) 
The real-space product in (4.1) becomes a convolution in Fourier space such that 



<^m2( k ) = 




where we used (3.21) to express the velocity perturbation in terms of the matter era amplitude 
of the Weyl potential, Ck- For large scales, k <C k s , it is easy to estimate the convolution 
using the approximation (3.20) for C^. The angular integration becomes trivial and the 
radial one can be split at k s . We find 

« *•> = l^(^r' (It) 2 irk + zk) < 44 » 

which provides a good estimate of the true numerical result that we compute using output 
from a modified version of the CMB-Boltzmann code CAMB [60, 61]. For the small-scale 
behaviour, k k s , we find analytically that 5p m2 drops off like ln(k/k s )(k/k s ) ns ~' i which also 
matches the behaviour of the numerical results. However, due to the involved angular integ- 
ration it is not straightforward to derive a simple expression for the amplitude on small scales 
and it is beyond the scope of this work to investigate this in more detail. In figure 1 we show 
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Figure 1. The matter pressure perturbation in the normal frame is shown for a typical ACDM 
cosmology at redshift z = (left panel) and z = 3 (right panel) as a function of k. The second 
order contribution from (solid red) falls below the leading third order contribution from & m \Vm\ 
(dashed black) around k « 0.4 h/Mpc at z = and around k ss 1 ft,/Mpc at z = 3. 
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Figure 2. The matter anisotropic stress in the normal frame is shown for a typical ACDM cosmology 
at redshift z — (left panel) and z — 3 (right panel) as a function of k. On large scales, k <C the 
term oc V^i (dashed black) is partly canceled by the contribution oc k~ 2 didj (V^V^) (dotted black) 
which is negative, such that their sum (solid red) falls off like k/k s . 



the numerical result for the effective matter pressure perturbation at redshifts and 3. To get 
an idea at which wavenumbers a second order calculation is not sufficient anymore, we also 
compute the leading third order correction, the term proportional to SmiV^. 4 The second 
order contribution falls below the leading third order contribution around k ~ 0.4 h/Mpc at 
z = and around k ~ 1 /i/Mpc at z = 3. 

For the anisotropic stress we take a double spatial divergence to construct a scalar 

4 We plot this third order term with the sole aim of showing, by comparison with the complete second 
order, at which intermediate (mildly nonlinear) scales, higher order correction terms can start to be relevant. 
However, our analysis is not meant to be a full perturbative treatment at those scales, which would probably 
require the use of resummation techniques (e.g. [62] and [63]). 
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Figure 3. The difference between the comoving and the normai frame matter power spectrum at 
second order in the longitudinal gauge, P* t — P m , is shown for a typical ACDM cosmology at rcdshift 
z = as a function of k. In the left panel the different contributions are plotted: (2u — 2/ 2 )P„ 
(solid red), 2P cross (dot-dashed magenta), v 2 — P au to (dashed blue). For comparison we show the 
linear matter power spectrum P^ 1 in dotted black. The right panel shows the relative difference, 
(P r *„ - P m )/P x ™, scaled by a factor of 10 6 . 



quantity 



Y2 d ^ 6 ^ =2 ^fJ (^fs [9 2 + 2k.q-3(k.q) 2 ]c k _ q C q (4.5) 



where k = k/fc. We compute the matter anisotropic stress numerically and show the res- 
ulting behaviour in k at redshifts and 3 in figure 2. When computing this convolution we 
realise that on large scales, k <C k s , the contributions from the term oc didj{y^i\Yml) * n (4-2) 
is negative and partly cancels the term oc V 2 ^^. We analysed the small- /c behaviour ana- 
lytically using a number of approximation schemes and found k~ 2 didjSTT l l ^ 2 (k <C k s ) oc k/k s , 
which correctly reproduces the numerical results. Again, an exact analytical estimate of the 
small-A; limit is difficult to achieve and not helpful for the discussion presented here. On small 
scales, k 3> k s , the behaviour of the anisotropic stress is dominated by the term oc V 2 !^ 
and therefore the scaling is also oc ln(k / k s )(k / k s ) ns ~ s as discussed above. 

Finally, let us study the difference between the density perturbations in the normal 
frame and the comoving frame. As we saw in section 2.2, the difference between the two 
point correlation functions of the density fields is due to the relative difference of the mean 
densities (2.9), and the pressure-pressure and pressure-density correlations, see (2.10). We 
compute this difference up to second order and express it in terms of the power spectra to 
find 



where 



P* m -P m = (2u - v 2 ) P l ™ - 2P CTOSS - [u 2 - P auto ] (4.6) 

(2vr) 3 Mk " *)F%(k) = <<W(k) Ci(k')) (4-7) 

(27rf5 D (k - k')Pcross(fc) = (<W(k) [V^VmuTQt')) (4.8) 

(2vr) 3 5 D (k - k')P au to(fc) = ([V^VmuKk) [V^V mll ]*(k')) . (4.9) 
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The linear densities coincide in the normal and comoving frames, and we find that v = 
Wmi^rnii\{k = 0) is the velocity dispersion. This gives a clear physical interpretation for the 
difference between the two frames: the geometric matter pressure seen by the normal frame 
is due to its velocity dispersion. In figure 3 we show the different terms in (4.6) and the 
relative difference between the matter power spectra in the two frames in a typical ACDM 
model. It is evident that the term 2i/P^ is the dominant contribution. In the present epoch 
we find the velocity dispersion to be v(z = 0) ~ 3 x 1CP 6 (for ACDM, in units c = 1, see 
figure 2), a value that corresponds to velocities of v « 500km/s. The fact that v is non-zero 
means that there is a fundamental difference between the comoving matter frame and non- 
comoving frames. Indeed, the mean densities are related by p* m = (1 — v)p m which means 
that the background evolution is slightly different in the two frames; different only at the 
level of second order perturbations. This means that in principle the background needs to 
be renormalised in order to take this into account, a procedure that belongs into the realm 
of backreaction calculations and is not of interest for the present work; see e.g. [64]. 

We conclude that the amplitude of 5p m today is about 10~ 6 p m on large to intermediate 
scales and is therefore small but a priori not negligible. We will compare it with the pressure 
perturbation in typical dynamical DE models in section 5.1. It is interesting to note that 5p m 
does not vanish for k — >■ which implies that, in principle, the large-scale average matter 
pressure is non-zero and needs to be renormalised. This was also recently pointed out in 
[64]. The matter anisotropic stress, on the other hand, is found to be vanishing for k — > 0. 
It peaks roughly at k ~ k s with an amplitude of ~ 10 p m . In the next section, we will 
compare it with the other contributions to the second order gravitational slip, II2 = </> 2 — 1P2, 
which form the second order consistency relations. 

4.2 Consistency relations at second order 

We calculate the geometrical fluid variables, 5p G2 , Sp G2 , 5q G2 and Stt 1 ^, as measured in the 
normal frame up to second order in the metric perturbations. The full expressions for all 
these quantities are presented in appendix C.3. Here we focus on the geometrical pressure 
and anisotropic stress, given in (C.15) and (C.17), that are used to write the second order 
consistency relations (3.3). Here we assume Hi = 0, but it is straightforward to include a 
non-zero first order gravitational slip. The anisotropic stress consistency reads 

97/ 2 

V 4 n 2 = — a 2 didjbTe G2 - 8^1 VV - 2(VV) 2 - 12V 2 (V9?i) 2 + 10(<9;d^i) 2 . (4.10) 

The second order gravitational slip, comes in as a source term in the second order Bardeen 
equation, which follows from the pressure consistency 

<f%,M o,„ m 3«5p G2 ,/V 2 H 2 Il>> H' 2 3 \ 

VV 7(V^i) 2 9 id (d) 2 , 

The geometric pressure perturbation and anisotropic stress are given by the effective matter 
stresses that we computed in the last section and any stresses coming from X^ u : 

SPG2 = S Pm2 + SPX2 , ^G2 = 67r m2 + 67T X2 ■ ( 4 - 12 ) 
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Figure 4. Different contributions to H/(2ip) are shown for a typical ACDM cosmology at redshift 
z = (left panel) and z = 3 (right panel) as a function of k. The 9 4 y 2 -type terms (dashed blue) 
mostly dominate over the term oc didj^Vmi^ml) (solid red), while on large scales, k < 0.01/i/Mpc 
the anisotropic stress from photons and massless neutrinos dominates (dot-dashed black). 

The consistency of ACDM requires 5p X 2 = = <5tt^ 2 . A few comments are in order: first, 
notice that the second order gravitational slip II2 is purely given by gravitational nonlinear- 
ities, C/ 4 c/? 2 -type terms, and the anisotropic stress measured by the normal frame, didj5it l Q 2 - 
This means that a gravitational slip is generated at second order regardless of the model, 
even in ACDM where there is no anisotropic stress in the comoving frame. Second, an ana- 
logous observation can be made about the source terms of the second order Bardeen equation 
(4.11). Again, apart from the gravitational nonlinearity terms, the source is given by the 
pressure and anisotropic stress measured in the normal frame. Finally, it is important to 
stress that a gravitational slip is generated at second order independently of our choice to 
work with the normal frame. We chose to use the normal frame because it provides a very 
direct way of arriving at the above equations and enables us to interpret the contribution of 
the stress-energy content to the nonlinear source terms. 

In the case of ACDM, we compute II2 (A;) from (4.10) by numerically convolving the 
potential computed with CAMB. In figure 4 we show the ^-dependence of the different 
contributions to 11/ (2ip) for a typical ACDM cosmology at redshifts and 3 where we compare 
it to the first order gravitational slip from photon and massless neutrino anisotropic stress 
(that is not very accurately computed in CAMB at late times and on small scales). The 
contribution from the <9 4 <^ 2 -type terms in (4.10) turns out to be larger than the contribution 
from the matter anisotropic stress, but comes with the opposite sign. Combining the terms, 
we find that we have to compute 



-1 



(27^ 



8kq 6 n + 5fcV(l + M ) ~ 6& W 



C|k-q|C<?- 



(4.13) 



For small k <C k s we find that the angular integration cancels the k and k s terms while the 
q integration results in a vanishing k 2 term. The leading order term picks up a logarithmic 
correction, and we find that <9 4 </? 2 oc k 4 ln(k/k s ), i.e. these terms go to zero somewhat faster 
than the matter anisotropic stress. 
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Figure 5. The three different relevant source terms of the second order Bardeen equation are 
compared for a typical ACDM cosmology at redshift z = (left panel) and z — 3 (right panel). The 
matter pressure source S$ p (dashed red) and the gravitational slip source Sjj (dot-dashed blue) are 
both positive, while the d 2 tp 2 -type source Sqi^i (dotted magenta) is negative. This leads to some 
cancellations in the total source term (black) which is positive (solid) on large scales, k < 0.01 h/Mpc 
(depending on redshift), and negative (long dashed) on small scales. 



Next, let us have a look at the second order Bardeen equation. The source term on the 
right hand side is a combination of terms from the matter pressure, the gravitational slip that 
includes the matter anisotropic stress, and further d 2 (p 2 -type terms. At weakly nonlinear 
scales the terms with time derivatives and those without derivatives can safely be neglected 
in comparison to the gradient terms. We write the total second order source term as 



S 2 
Sn 

'<9V 



Ssp + Su + Sqi^i 

v 2 n 2 
m 2 



ml 

1 + cr 3 



30n(l + a 3 



n, 



VVi , 7(V^i) 



m 2 m 2 

aJ-(-a 3 ) 2 2 
2fUl + a 3 ) 



d 3 q 
(2vr) 3 



((? 2 + 7k-q)C k _ q C q . 



(4.14) 
(4.15) 

(4.16) 

(4.17) 
(4.18) 



We compute the convolution in Sq2^2 numerically using <p(k) from CAMB for a typical 
ACDM cosmology at redshifts and 3 and compare it to the other two source terms in figure 
5. We find that 5 , g2 y ,2 is negative while S$ p and Su are positive, such that some partial 
cancellations occur in the total source term, S2. This turns out to be positive and roughly 
constant on large scales, k < k s , where it turns negative and starts to decay in the same 
way as Sp m2 , so oc ln(k/ k s )(k / k s ) ns ~ 3 . Using the source term derived here, it is now in 
principle possible to solve the second order Bardeen equation for 932, but we leave this for 
future work since it is not the aim of this paper. We emphasise however, that the derivation 
of the Bardeen equation was simple taking the nonlinear consistency relations as a starting 
point. It is straightforward to generalise this to higher order. Moreover, this makes clear 
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how a possible dark energy pressure, 8p X 2i would come in as a source term. Due to the 
subtle cancellation we find, it is likely that the dark energy pressure would dominate the 
source term on small scales if the dark energy sound speed is not too close to unity and 
its equation of state parameter is not too close to —1. If this is the case, the sign-change 
and drop-off in S2 moves to smaller scales, see figure 5, and (f2 grows instead of decaying 
where the source becomes positive and appreciable. Therefore we can expect tp2 to strongly 
depend on the dark energy model. Extrapolating from linear theory, we expect that specific 
combinations of the Weyl potential, cp2, and the gravitational slip, II2, are directly linked to 
WL, peculiar velocity and galaxy clustering observables via the Euler and Poisson equations. 
Therefore, their behaviour on mildly nonlinear scales can provide an interesting smoking gun 
for non trivial dark energy dynamics that to be tested by combining different cosmological 
probes. A determination of the dependence of cp2 on the dark energy properties and, hence, 
of its effects on the aforementioned observables, cannot generally be done without resorting 
to specific models, but the argument above nevertheless shows that these may be used as 
probes of dark energy. 

5 Applications and discussion 

5.1 Dark energy pressure perturbation 

In ACDM the matter pressure 5p m from (2.5) is the only source of pressure perturbations well 
after radiation-matter equality. However, if the late time accelerated expansion is not driven 
by a cosmological constant but by some form of dark energy in GR, we can expect an intrinsic 
pressure perturbation 5p x to appear along with fluctuations in the dark energy component. 
In this case, according to the first equation in (3.3), the total pressure perturbation Sp G 
is the sum of 5p m and Sp x . This raises the question of which one of these two terms is 
more important for generic dark energy models. Whereas the first non-zero contribution to 
matter perturbations (3.13) appears at second order, dark energy perturbations are already 
present at first order. However, dark energy density perturbations are typically expected 
to be small in comparison to matter ones, S x ~ (1 + w x )5 m , and the behaviour of the 
dark energy pressure perturbation depends on the sound speed. Therefore, it is a priori 
unclear how matter and dark energy pressure perturbations compare to each other. Here, 
we compare the pressure perturbation of a perfect fluid dark energy model with the effective 
matter pressure that we quantified in the last section. We do this to find out whether the 
matter pressure could, in principle, be mistaken for a dark energy pressure perturbation, or 
if it has to be taken into account when trying to reconstruct the dark energy properties from 
its pressure perturbation. As we will see next, the precise answer to this question depends 
on the wavelength of the fluctuations. 

It is well known that the amplitude of dark energy fluctuations increases as the equation 
of state parameter w x deviates from —1. In addition, the importance of the dark energy 
clustering is determined by the sound speed which sets the Jeans length ~ c x Tl~ l below 
which pressure support inhibits gravitational collapse. In the limit in which the sound speed 
equals the speed of light, clustering becomes negligible and therefore virtually impossible to 
detect. For an incomplete list of references where dark energy perturbations are studied see 
[65-76]. 

The sound speed of dark energy can be defined as the relation between its pressure and 
density perturbations in the rest frame (corresponding to the comoving orthogonal gauge), 
5p x = c x 2 5p x . It captures possible non-adiabatic pressure perturbations and therefore 
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Figure 6. Comparison of the dimensionless power spectra of dark energy pressure perturbations 
and effective matter pressure Sp x /(p m + px) (black & blue) and 8p m /{p m + px) (red), respectively, 
at z = (left panel) and z — 3 (right panel) . For the matter pressure we show the second order result 
oc V^i (solid red). For the clustering DE we compare two models: Wx = -0.8, c x 2 = 0.001 (solid 
blue) and wx = —0.98, c'x 2 = 0.01 (dashed black). The equation of state, wx, controls the amplitude 
while the sound speed, c3f 2 , controls the scale at which the spectrum starts falling off on small scales. 
The dotted asymptotes to the dark energy spectra for small and large scales are given by the analytical 
approximations (5.3) and (5.4). The intrinsic dark energy pressure is orders of magnitude larger than 
the effective matter pressure for large scales. However, for small scales, k > 0.3/i/Mpc, the second 
order matter pressure starts to become relevant depending on the dark energy model parameters and 
rcdshift. 



is, in general, not equal to the ratio p'x/p'x- The rest-frame sound speed, cx, defines a 
relation between gauge-invariant combinations that can easily be transformed into any gauge 
[6, 7, 65, 77]. Assuming for simplicity a constant equation of state we have 

$Px = cx 2 5p x + 3(1 + w x )(c x 2 - wx)-^PxVx (5.1) 

which is valid in any orthogonal gauge (and for any frame). We will focus on the case of small 
(positive) values of c'x 2 , which is the most interesting case observationally. Then, the first 
term on the right hand side of the expression above is negligible on large scales, k <C H/cx, 
and the dark energy component is approximately comoving with dark matter, Vx — V m , see 
e.g. [72]. So we can write: 

8px-Swx(L + wx)^rPxV m , iork<^H/c x (5.2) 
where the dark matter velocity can be read off from (3.21) and (3.22). Using the properties 
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of the hypergeometric function 2F1, this result can be approximated at z = by 

6 Px (z = 0)~-^w x (l + w x )p x (l + -^-\c k , for k<H/cx (5.3) 

where Ck for k <C k s is given in equation (3.20) in terms of the primordial power spectrum. 

On small scales, deep inside the nonlinear regime (and hence below the sound horizon 
of dark energy), the relevant term of (5.1) is the one that depends on Sp x . In this limit it 
turns out that 5p x is inversely proportional to c x 2 [69]. Therefore the pressure fluctuation 
is again nearly independent of the sound speed 

$Px - -{l + w x )p X (p~ -{l + w x )p x fe- ^a~ 3 ^) a^Ck, for k^>H/c x (5.4) 

where the Weyl potential ip = C k F{—a^) has been approximated by its late time expression 
(C.7) (where also e « 1.44 is defined), while Ck is given by the k 3> k s branch of (3.20). 
These results show that the amplitude of the linear dark energy pressure perturbation does 
not depend on the value of the sound speed either for large or small scales. Naively, one 
would expect that nonlinear corrections would modify this feature on small scales. However, 
the results of [74] on the density power spectra indicate that the nonlinear corrections are 
roughly independent of the (small) value of the speed of sound. 

What the sound speed does affect is the behaviour of the pressure perturbations at 
intermediate scales. This is shown in figure 6 where we plot the dimensionless power spectra, 
as defined in (A. 3), of Sp m /(p m + p x ) and Sp x /(p m + p x ) at redshifts and 3. Notice that 
we divide the pressure perturbation in each case by the total energy in the clustering DE 
model. Here we assume that 5p m computed in ACDM is roughly equal to 5p m in clustering 
wCDM. This is well justified as 5p m only depends on V m which is not affected strongly by 
dark energy clustering. The two models of dark energy depicted assume constant values of 
w x and c x 2 . The amplitude at small and large scales is given by w x and is independent of 
c x 2 . Hence, the equation of state controls the overall amplitude of the pressure perturbation 
and the sound speed determines the scale at which 5p x starts to drop off. 

The power spectrum of 5p m is computed at second order in the velocity perturbation. 
On large scales, the dark energy pressure perturbation is much larger than its counterpart 
from the frame effect on the matter. However, even though a complete nonlinear analysis is 
needed to study small scales (k > 0.1/t/Mpc) the results indicate that for nonlinear scales 
the matter pressure can be comparable to the dark energy one. Therefore, for studies of 
dark energy perturbations on large scales using, for instance, the ISW effect on the CMB, 
the effective matter pressure can be safely ignored for interesting values of the dark energy 
parameters. On the other hand, for studies of dark energy clustering on smaller scales (e.g. 
WL) , the frame effect will have to be taken into account for a reliable interpretation of the 
data. Such a study is beyond the scope of the present work and will involve an accurate 
computation of the dark energy perturbations on nonlinear scales [74, 75] including also the 
higher order contributions to (4.1). 

5.2 Anisotropic stress in modified gravity 

While the existence of a sound horizon and the associated pressure perturbations are a key 
feature of scalar field dark energy models like Quintessence and K-essence, the anisotropic 
stress is usually considered as a smoking gun for modifications of gravity [67, 78, 79]. However, 
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we have shown above that at nonlinear scales even ACDM exhibits a non-zero gravitational 
slip. This raises the following question: would it be possible to confuse this higher-order 
contribution with the anisotropic stress intrinsically present at the linear level in some MG 
models? Here we consider the DGP model which assumes that we live on a 4D brane 
embedded in a 5D bulk. Even though the self-accelerating solution has a ghost degree of 
freedom [80, 81] and, in addition, this model is observationally ruled out at the 5a level (e.g. 
[82, 83]), it is nonetheless one of the 'standard' MG models as it exhibits many of the typical 
features of this class, and it has been well studied in the literature. Hence, for our purpose 
of comparing the matter anisotropic stress due to the frame effect to the anisotropic stress 
in MG, it will serve us well. 

The gravitational DGP action is just the 5D Einstein-Hilbert action together with an 
induced 4D Einstein-Hilbert action confined to the brane. The relative strength of the two 
contributions is given by the crossover scale r c = Mp/Mf where M5 is the 5D Planck mass. 
DGP therefore, like ACDM, is characterised by a single free parameter, the crossover radius 
r c . The background expansion (assuming flatness) in DGP is given by (see e.g. [84]) 



R 2_H = Jrn ( } 

r c 3MJ, V 7 



where H is the cosmic time Hubble parameter. We will consider the second term on the 
left, H/r c , as providing the dark energy contribution to the expansion rate and so think of 
it as being effectively on the right hand side. If the Hubble parameter today is given by 
T-Lq and the relative matter density today is taken to be Q m , then we want for consistency 
7~io/ r c = HqO- — ^m) or r c = 1/[Hq(1 — f2 m )]. As naively expected, we see that the crossover 
scale of DGP needs to be of the size of the horizon scale today, which reveals that the 
fine-tuning required in this model is comparable to the one in ACDM. In the quasi-static 
limit relevant for scales that are subhorizon but still linear (see e.g. [85]), one finds for the 
perturbations that 

- k ^= a ^X i -^^ + ^ Vm ) (5 - 6) 

= f * + i V<W + ^V m ) (5.7) 



2Mp V ' 3{3 J ' k 

for /3 = 1 — 2r c H[l + H / (3H 2 )] = 1 + 2r c Hwx, and the matter perturbations evolve as usual 
as a pressure-free fluid sourced by the gravitational potentials above. However, since the 
Poisson equations for the potentials are altered, also the matter evolution is now different. 
Defining 771, the MG parameter related to the anisotropic stress, as in [26] we find that 

III 2 , , 

'^-urrwi- <5 ' 8) 

Numerically, rji is small at high redshifts but tends towards rji ~ —0.44 today for a flat Uni- 
verse with Q m = 0.3. In other words, the effective anisotropic stress in DGP on linear scales 
is of comparable size as the gravitational potentials themselves, and much bigger than the 
second order anisotropic stress in ACDM. Current constraints on the MG parameters are very 
weak, with deviations of order unity still allowed even for fairly restrictive parameterisations 
(see e.g. [12, 32, 33, 35] for recent analyses). 
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Once the matter perturbations become nonlinear, the situation is more complicated. 
All viable MG models, specifically all models that pass solar system tests, need to screen the 
extra degrees of freedom at small scales since they would induce a fifth force of a strength 
far beyond observational limits. This screening effect is usually coupled to the local density 
through a Vainshtein [86], Chameleon [87] or similar mechanism (e.g. [88]) and thus becomes 
active on small scales. In the case of DGP it is the Vainshtein mechanism that provides the 
necessary screening, and it is generically active on scales of the order of [89, 90] 



where r$ = 2GM is the Schwarzschild radius of the source (here assumed to be a point mass). 
For a cluster with M = 1O 15 M0 the Vainshtein radius is of the order of r* ~ 10 Mpc, which 
indeed coincides roughly with the scales at which perturbations start to become nonlinear. 
On smaller scales, ref. [91] estimates the anisotropic stress to be 



which means it is strongly suppressed on small scales, and roughly goes as —(r c k)~ 2 (f)i. 

To summarise the discussion, the power spectrum of IT in the DGP scenario has the 
following order of magnitude today: 



The anisotropic stress is therefore large on scales above roughly 0.1/Mpc and decays more 
rapidly on smaller scales. The exact behaviour would have to be studied with numerical 
simulations. Comparing the small size of the matter anisotropic stress relative to the Weyl 
potential shown in figure 4 to the order unity effect in DGP shows that on linear and weakly 
nonlinear scales the MG signal is not contaminated by the matter velocity induced gravita- 
tional slip. The difference is so large that we expect the dynamical anisotropic stress in MG 
models to remain smaller than the effective matter anisotropic stress on all scales of interest 
for future large-scale surveys like e.g. Euclid 5 [92]. 

5.3 Confusion with galaxy bias? 

Finally, there is the question of whether it is possible to construct combinations of observables 
that are sensitive to the effective anisotropic stress and pressure. And, if the answer is 
positive, can we disentangle those signals from other uncertainties? A good example is 
the reconstruction of the galaxy bias by cross-correlating galaxy number counts and WL 
observables [41, 93, 94]. 

For illustration, let us briefly discuss this in linear theory with a slight gravitational 
slip, II = (j) — ip 0. Our aim here is to clarify how the presence of a gravitational slip, if 
ignored, would be misinterpreted as a non-trivial galaxy and velocity bias. Well below the 
Hubble scale, the observed galaxy number density fluctuation in the direction n is 




(5.9) 




(5.10) 




k < 1/r* 
k > 1/r* . 





(5.12) 



5 Euclid consortium webpage: http:// 



www . euclid-ec . org 
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in terms of the actual galaxy number density fluctuation in the longitudinal gauge, 5 g , the 
galaxy peculiar velocities, V g , and their WL convergence, k (see for instance [95-98]) 6 . Here 
Hk = rik. We apply a simple local model for the Lagrangian galaxy bias where the number 
density, g, is a function of the comoving (rest frame) matter density, p^, and the proper 
time, r. In the longitudinal gauge (normal frame) we can then write 

5 g = b5 m + 3bjV m + jV g , V g = b v V m (5.13) 
where we define the Lagrangian galaxy bias b and the evolution bias b T 

as derivatives of the average galaxy number density g. The velocity bias b v is often taken to 
be unity but we keep it here for clarity. 

We now write all the perturbations in terms of the metric potentials. On subhorizon 
scales, and neglecting time derivatives of the potentials, we have from (3.14) and (3.16) that 

(* + i n ) (5 ' 15) 

The WL convergence of sources at conformal distance Xs is 

rxs 

«(n,Xa) = V|/ dxW{x s ,x) l p(Xsn,r] - x) (5.17) 
J o 

where V| is the Laplacian on the unit sphere and 770 is the conformal time today. The WL 
efficiency function in a flat background is W(xsiX) = (Xs — x)/(XsX)- Its effect is a very 
broad 3D to 2D projection while weighing down contributions of lenses close to the sources 
as compared to those close to the observer. Therefore it is convenient to cross-correlate 
the galaxy number counts at a given redshift, A g = A(n,z g ), with the WL signal of the 
background sources at higher redshift, Kb = «(n, Zb), such that the galaxies select a thin slice 
from the broad WL efficiency function. Well inside the horizon we can safely neglect the 
velocity term otH/k in (5.13) such that 

(A 5 Kb) ^ {b5 m K b ) - T-r l {b v p k kV m K b ) - 2{nK b ) . (5.18) 

Using the expressions for 5 m and V m above, we can write schematically 

(A g k„) (k 2 [b(i + n/2g) + b v (i - n/2y)] g y g w m ) 

where the line-of-sight integrations of the WL convergence are implicitly assumed. It is now 
easy to see that the presence of anisotropic stress, or a gravitational slip, leads to replacing 
the standard biases by 

6 - 6 (l + °-), b v ^b v (l-^). (5.20) 



Here we neglect the corrections from the Doppler term, the ISW effect, the time delay and the metric 
potentials at the galaxy positions for simplicity and also because they are negligible well inside the Hubble 
horizon. 
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Or the other way round: ignoring the presence of a gravitational slip biases the measurement 
of galaxy biases as shown above. Using different tracers of the density field at different 
redshifts one can now disentangle the biases and the gravitational slip. The corrections from 
second order perturbations are very small, see figure 4. However, this result holds for any 
metric theory of gravity and most MG models do generate significant anisotropic stress as 
discussed in the previous section. Therefore it is in general important to take this effect on 
the bias into account. 

In [41], an elegant non-parametric approach was developed to test GR on cosmological 
scales using a set of consistency relations involving metric and matter density perturbations 
in linear theory. It was pointed out that to test any of these relations, the observations 
of peculiar velocities and galaxy clustering can be combined to predict the outcome of WL 
measurements in order to have a strong self-consistency test of ACDM and measure the linear 
galaxy bias similarly as discussed here. While this proposal is consistent at linear order, we 
would like to point out that, in any case, nonlinearities induce anisotropic stress that needs 
to be taken into account properly when undertaking such a program with information on 
nonlinear scales. Furthermore, as we see above, if a gravitational slip is actually present in 
linear theory in reality (as in MG), it is readily mistaken for a non-trivial galaxy bias if not 
taken into account. 

6 Conclusions 

With current cosmological observations tightening constraints on the ACDM scenario and 
not (yet) hinting at strong deviations from it, or from GR, it becomes increasingly important 
to devise ways of falsifying or confirming the standard paradigm. To be able to do so, a 
very precise understanding of the observations and the corresponding predictions is needed 
not only on linear but also on nonlinear scales. In recent years, much progress has been 
made in nonlinear Newtonian perturbation theory to understand and complement numerical 
analyses. Simulation efforts in non-standard scenarios are also underway. On the other hand, 
relativistic perturbation theory predictions and consistency relations are mostly restricted to 
linear theory. But to be able to confidently falsify or confirm the standard paradigm it is 
imperative to understand the implications of dynamics in curved space on weakly nonlinear 
scales where Newtonian theory may not be sufficient, even more so if modifications to GR 
are allowed. 

A framework based on casting beyond-ACDM structure growth in terms of two func- 
tions, that can either be parameterised in generic ways or predicted from a given theory, 
has recently received considerable attention. This tool is often referred to as Parameterised- 
Post-Friedmann approach or modified growth parameterisations. In this paper, we propose 
an alternative framework deriving fully nonlinear consistency relations between the geometry, 
the matter content and any possible dark component. Since our framework requires neither 
the assumption of GR nor minimal interactions, but only relies on the symmetries of the 
Einstein tensor, it can be applied to different scenarios and models. The framework derives 
from the covariant approach, employing a very practical choice of 4-vector field (the normal 
frame) to decompose the Einstein tensor, the energy-momentum tensor and their difference, 
X^ u = M P G^ V — T^ v ', that describes the physics beyond GR and/or beyond the standard 
constituents of the cosmological model. The normal frame, n^, is defined to be orthogonal 
to the surfaces of constant time. Because its definition is purely geometric it is only given by 
the metric itself and so are the irreducible components of MpG^ u : the geometrical dynamic 
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quantities. This is the reason why the normal frame is such a convenient choice, as it natur- 
ally allows us to write consistency relations that separate into a geometrical part, a matter 
part and a dark energy part: 

1 2 ixv iiv (u u) 

PX=PG~ ^mPm , =7T G ~ Pm v m v m 

where the G- variables are only given by the metric, see section 3.1. These fully nonlinear 
consistency relations can either be taken as conditions for ACDM+GR to hold by requiring 
px = —MpA, n't 11 = 0, or the X-quantities can be parameterised and reconstructed from 
data. The advantage of this very general framework is that it is a top-down approach that 
a priori takes all non-linearities and all scalar, vector and tensor degrees of freedom into 
account, and can be simplified in perturbation theory at will. 

As a consequence of selecting the normal frame, the dynamical quantities on the matter 
side (the energy density, the pressure, the energy flux and the anisotropic stress) do not take 
the values that one would expect, or find when measured in the matter rest frame. The 
transformations from the rest frame to the normal frame can be written fully nonlinearly. 
These show that a non-vanishing effective pressure and anisotropic stress are measured in 
the normal frame, even if the matter is assumed to be a perfect pressureless fluid in its rest 
frame, see section 2.2. This is the origin of the terms ~ in the consistency relations. 
However, these effects are only relevant at second order in perturbation theory and they can 
be accounted for, knowing that they are present. Furthermore, we point out that cosmolo- 
gical observations do not actually measure the density field in its comoving frame, as we do 
observe the effect of the velocity field of the sources through the projection onto the photon 
wave vector. This means that there are non-comoving effects that need to be taken into 
account. More importantly, apart from gravitational nonlinearities, we identify the pressure 
and anisotropic stress measured in the normal frame as the physical quantities that enter as 
sources in the Bardeen equation for the second order Weyl potential and the constraint for 
the second order gravitational slip. This is a consequence of the fact that the normal frame 
is purely defined in terms of the metric. Finally, let us recall that the normal frame helps to 
recover the Eulerian picture of Newtonian perturbation theory from GR. Although GR is a 
theory of curved spacetime, the normal frame is related to the Eulerian picture in Newtonian 
nonlinear dynamics. 

We compute the consistency relations for ACDM in first and second order scalar per- 
turbation theory in the generalised longitudinal gauge (that is the gauge where n^ 1 has zero 
spatial components and does not see shear or vorticity). At first order we recover the usual 
conditions: the gravitational slip vanishes, H\ = (pi — ipi = 0, and the Weyl potential, 
ipi = \{<p\ + ipi), follows the standard Bardeen equation, see section 3.4. We solve the 
Bardeen equation analytically to provide the basis for a careful study of the matter pressure 
and anisotropic stress. We compute ^v^Pm arid p m VmVm at second order and find these 
terms to be small on large scales. They become relevant only on scales k > 0.1 /i/Mpc. The 
second order anisotropic stress consistency relation reveals a non-vanishing gravitational slip 
due to the effective matter anisotropic stress and non trivial combinations of products of spa- 
tial derivatives of (p±, see equation (4.10) and also refs. [99, 100]. Here we clarify the origin 
of the different contributions to Tl2- The pressure consistency, on the other hand, yields the 
second order analogue of the Bardeen equation for ip2, see equation (4.11). It acquires differ- 
ent source terms which we analyse, compute and compare and we find subtle cancellations 
between the contributions from the effective matter stresses and d 2 (f1-type terms. The fact 
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that the source term of the second order Bardeen equation depends non trivially and sens- 
itively on the pressure and anisotropic stress perturbations of a possible dark energy makes 
the weakly non-linear evolution a promising tool to test the ACDM paradigm. Combining 
different cosmological probes enables us to reconstruct the consistency relations, as all ob- 
servations probe the metric potentials, either directly or indirectly via the Poisson equation, 
see also e.g. [41]. Let us emphasise that the effective matter pressure and anisotropic stress 
are seen by the second order metric perturbations in any case, irrespectively of the yet to be 
solved issue of which precise frame is the one that corresponds to us as observers. 

Finally, we discuss a set of applications: we compare the effective matter pressure with 
the pressure perturbations in clustering dark energy, concluding that the matter pressure is 
negligible on large scales and becomes relevant on weakly nonlinear scales. Thus, the matter 
pressure only needs to be taken into account when deriving constraints on the dark energy 
sound speed from small scale observations. Furthermore, we compare the matter anisotropic 
stress to typical anisotropic stresses arising in MG scenarios, specifically the DGP model. 
Anisotropic stress is normally considered a smoking gun for the detection of MG. We find 
that, typically, the matter anisotropic stress is by far subdominant and not an issue. Finally, 
we assess the impact of the anisotropic stress on the determination of galaxy bias from 
combining WL and galaxy clustering data. We find that the gravitational slip affects the 
Lagrangian galaxy bias and the velocity bias in opposite ways, but only marginally. Using 
multiple tracers at several redshifts and different combinations of observables, we expect that 
one could be able to filter out or strongly constrain the gravitational slip. 
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A Conventions 

We employ natural units where c = h = kp = 1 and define the reduced Planck mass as Mp = 
(87rG) _1//2 , with G being Newton's constant as measured in laboratory experiments. We 

work with a signature ( h + +)• For tensor components, Greek indices take values ... 3, 

while Latin indices run from 1 to 3. Round parenthesis around indices mean symmetrisation, 
y(H = i(y^ + y^) ; an d squa re brackets mean anti-symmetrisation, = \{Y liU -Y u ^). 

The spatially flat, homogeneous and isotropic FLRW metric is g^ v dx^dx u = a 2 (—dr] 2 + 
5ijdx l dx J ), where r/ is the conformal time and a{rf) is the scale factor normalised to unity 
today so that the comoving scales are numerically equal to physical scales, at present time. 
Derivatives with respect to rj are denoted with a prime and the conformal Hubble parameter 
is T~L = a! I a = aH where H is the cosmic time Hubble parameter, whose value today is the 
Hubble constant Hq = 100/ikm/s/Mpc rj 3 x 10~ 4 /i/Mpc. Overbars denote background 
quantities and subindices in perturbation variables indicate the order of the fluctuation. In 
our notation, represents the present relative energy density of the i-component, and is 
therefore a constant unless specified otherwise. 
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When numerically computing results for a typical flat ACDM model we use the following 
parameter values: A s = 2.1 x 1CT 9 , n s = 0.96, k p = 0.05 /Mpc, h = 0.7, Q c h 2 = 0.112, 
tt b h 2 = 0.0226 and fl r h 2 = 4.17 x 10" 5 . 



bfi = u.uzzd ana ii r t 

Our convention for the Fourier transform is 



/(k) = J d 3 x /(x)e lk ' x , and hence /(x) = J /(k) e"** . (A.l) 

The power spectrum of any random variable / is given by the Fourier transform of the 
two-point correlation function: 

(2vr) 3 Mk " *)Pf{k) = (/(k)/* (k')) (A.2) 
where <5o(k — k') denotes the 3D Dirac delta distribution. Finally, 

Vf{k) = ^Pf(k) (A.3) 
defines the dimensionless power spectrum. 

B Covariant 1+3 decomposition 

This appendix gives a brief summary and an account for our notation of the covariant ap- 
proach to cosmology. More details on the formalism and its applications can be found for 
example in [42-47]. 

B.l Defining the covariant formalism 

The covariant approach defines a 1+3 splitting of space-time by using a time- like unit 4- 
velocity field, with u 2 = g^ u u^u u = — 1. It represents a family of observers and is often 
referred to as an observer frame. All tensors are decomposed into their irreducibles w.r.t. u M , 
the projection tensor = g^ v + u^u v and the projected alternating tensor e^ ua = n^apu 13 
where ?/ M!/Q( g = — \J— g89 S^d^S^ . Then a dot and represent the covariant derivative 
parallel and orthogonal to u^, respectively, which for a generic tensor J read 

>•"...„ = »"V„./""... ( , (B.l) 
D x J^... v = V%V$ ■ ■ ■ nV a J P '"..^ (B.2) 

To describe the geometry as seen by the family of observers u^, one uses the irreducible 
components of V^u^ 

= -ii^Uu + 3®?V + £ ^aUJ a + (Tfu, ■ (B.3) 

These are the kinematic or geometric quantities of u^: the expansion = D a u a , the accel- 
eration = u a V a u^, the shear a^ v = D^u^, and the vorticity u:^ = —\e^ a pD a u^. Here 
the projected symmetric trace-free parts are defined as 



Y aP . (B.4) 



In the covariant approach we decompose Y^ v 6 {MpC 1 ", T^, X^ u } into its dynamic 
quantities or fluid variables: energy density py = UaUpY"^ , isotropic pressure (including 
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possible bulk viscous stress) py = \V a pY a P , energy flux (heat plus particle flux) q Y 
—u a Y^ a and anisotropic stress (shear viscous stress) 7T Y U = Y^', Then, we can write 



Y" u = PY u»u v + pyV^ + 2q^u^ + tt y u 



(B.5) 



and the dynamic quantities add up as follows, see (3.1), 



PG = PT + Px , PG = Pr + Px , q G = q T + ( ixi 



G 



uv , uv 

ITrp + 7Ty 



(B.6) 



where the subscript G indicates the "geometrical" fluid variables defined by projecting 
MpG^ u . Notice that this is still fully nonlinear and generic, and these are dynamical and 
constraint equations for the metric and the possible additional dark degrees of freedom, but 
there are no derivatives of the matter variables coming in. Usually we see the G-variables 
simply as the total energy density, pressure, etc. 

B.2 Change of frame 

The fact that each component has an energy flux is because we decompose the multi- 
component system w.r.t. a generic observer frame rather than a specific 4-velocity of 
one of the components or the total energy- momentum content. However, the physical prop- 
erties such as the sound speed or the barotropic equation of state of a given component are 
described in its energy-frame, which means in terms of the dynamic quantities as observed 
by its 4-velocity u Y , denoted by an asterisk. The energy- frame is the frame in which the 
energy-flux vanishes, qy 1 = and in most cases coincides with the particle rest frame. The 
energy- frame satisfies the eigenvalue equation YaU Y = —p Y u Y . In terms of the energy- frame 
dynamic quantities the energy-momentum tensor reads 



Y^ V = PyllyUy + PyTy + 7Ty 



*flU 



(B.7) 



g^ v + UyUy. The transformation from the rest frame u Y to a generic frame 



where V Y V 

depends on the relative velocity of the two frames, v Y . It is given through the relation 
u y = 7y(u^ + v Y ) that is derived by decomposing u Y w.r.t. and using the normalisation 
conditions u Y = —1 = u 2 . Therefore v Y u^ = and jy = (1 — v Y )~ 1 / 2 . The dynamic 
quantities in the rest frame transform to the generic u^-frame as 



PY = Py + 



Py = Py + 



Qy 



TT' 



I-"' 
Y 



2 2 i * i * \ i *"P 
lY v Y (PY + PY) + ^Y v YaVY/3 



a/3 



2 2/* i * \ i *aB 
lY v Y\PY +PY)+ 7r Y v YaVYj3 



(Py +Py) v y+ 7y v y (Py +Py) v y + Ky^VYc 



*a3 a 
7T Y VYaVYfiU 



TX Y + 



1y(Py+Py)v^v y } 



Y V Y 
1 



(B.8) 
(B.9) 
(B.10) 

(B.ll) 



see for instance the appendix of ref. [47] specialised to u Y being the rest frame of the Y- 
component, i.e. q^ = 0. The terms in square brackets are at least second order in the relative 
velocity and only the energy flux is non trivial at first order. 



-27- 



B.3 Conservation equations 

The connection between kinematic and the dynamic quantities is made via V^G^ = 
a consequence of the Bianchi identities and therefore a geometrical fact (or total energy- 
momentum conservation if you wish): 

PG + e(p G + PG ) + (D a + 2u a ) q G + (r ai g7r£f = (B.12) 
€) + 3 ^ + MPG + VG) + D^po + (a m + e^J*) q G + (D a + u a ) Trg, = (B.13) 

which are the continuity and the Euler equations, respectively, and are not independent 
of the field equations above. Note that the total G-variables can only be replaced by the 
T- or X-variables under the assumption that the respective energy-momentum tensors are 
separately conserved. 



C Metric perturbations and the normal frame 

C.l The kinematic quantities of the normal frame in the longitudinal gauge 

The normal frame defines a set of kinematic quantities in the covariant approach: acceler- 
ation n M , the isotropic expansion 6, the shear a^ u and the vorticity u^. Up to second order 
in cosmological perturbations in the generalised longitudinal gauge they read 

1. 
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(C.l) 
(C.2) 



(C.3) 

(C.4) 
(C.5) 



Their covariant definitions are given in appendix B.l. 

C.2 Exact solutions of perturbation evolution after radiation-matter equality 

The Bardeen equation after radiation-matter equality (3.18) can be solved analytically in 
ACDM where w G = -0 A /(0 A + O m a -3 ). We find 

Vl = C7 k ^(_ a 3) + zWa~ 5 + a" 2 (C.6) 

where a = (f}A/^m) 1//3 « is the scale factor normalised at A-matter equality. The second 
mode, oc D^, is decaying rapidly and can safely be ignored. The time-evolution of the first 
mode, oc Ck, is given by the ordinary hypergeometric function F{y) = 2^(1/3, 1; 11/6; y). 
To see its early and late time behaviour we expand the solution for small and large a to see 
the limiting behaviour: 
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IT 
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l' 



-a 



1 



-of 



-a 3 ) a ] 



for a <C 1 
for a 3> 1 



(C.7) 
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where e = 27r~ 1 / 2 r(2/3)r(ll/6) ~ 1.44. The early-time approximation is accurate to < 1% 
for a < 0.5 and to < 10% for a < 0.8, while at a = 1 the late time approximation is more 
accurate than the early-time approximation, although only good to about 22%. 

We derive the matter density contrast and the velocity perturbation from the exact 
solution for (p\ by means of (3.14) and (3.16): 



v ml 



2a tr -2(l + a 3 ) 
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din a 



<Pl + 



d(pi 
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(C.8) 
(C.9) 



where we used pc/ Pm = (1 + m ACDM. Using the solution (C.6) with = we find 
for the matter density contrast 



Jml,k 



aF{-a 3 ) k 2 



3n r 



,, 2 + (l + a 3 ) T(-a 3 ) - 3a 3 T(-a 3 ) 



(CIO) 



where F(y) = {d/ dy)F{y) = 2-^i(4/3, 2; 17/6; y). On subhorizon scales, k 3> "H, we recover 
the usual growing mode solution, 5 m i oc a, in the matter era, using the early-time behaviour 
of T given in (C.7) for aCl. For the velocity, we find 
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which is also valid on all scales. At early times we have 
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(C.ll) 
(C.12) 

(C.13) 



an approximation that is accurate to < 1% for a < 0.4 and to < 10% for a < 0.6. 

C.3 Second order geometrical fluid variables after radiation-matter equality 

We compute the geometrical fluid variables, Sp G2 , <5pg2> $Qg2 anc ^ ( ^ 7r G2 a ^ secon d order in 
the generalised longitudinal gauge. These are the irreducible dynamic variables projected 
from MpGpv, see section (B.l). To simplify the expressions, we already impose the ACDM 
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condition for the first order gravitational slip, 111 = 0. We find 
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These expressions are used in section 4.2 to derive the consistency relations at second order. 
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